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Abstract. Vacancies in graphene lead to the appearance of localized electronic states 
with non-vanishing spin moments. Using a mean-field Hubbard model and an effective 
double-quantum dot description we investigate the influence of strain on localization 
and magnetic properties of the vacancy-induced states in semiconducting armchair 
nanoribbons. We find that the exchange splitting of a single vacancy and the singlet- 
triplet splitting for two vacancies can be widely tuned by applying uniaxial strain, 
which is crucial for spintronic applications. 
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1. Introduction 

Spintronics is an interesting prospective technology where spin is nsed as the information 
carrier (in contrast to electronic charge in the case of electronics) jU [2]. Graphene has 
been proposed as a promising material for realizing this technology |3]. A small spin- 
orbit conpling facilitates long spin decoherence times and makes information processing 
based on spin feasible [H |5]. Magnetic order can can arise in graphene nanoflakes[6l 0 [HI 
and in connection with localized edge states in zig-zag nanoribbons|inj. Farther, point 
defects in graphene, snch as vacancies or chemisorpted hydrogen, have been proposed 
for realizations of spintronic devices fTTlfT^] . since those defects snpport (qnasi) localized 
electronic states and can indnce magnetism [laiiiiiisi IT6] . In practice, snch defects 
can be created by nsing electron or ion beams fl^ FIT] . 

The magnetic order associated with the magnetic moments depends on the spatial 
conhgnration of the defects nam]. For gapped graphene strnctnres, snch as nanoflakes 
[H] or armchair graphene nanoribbons (AGNRs) [TH], the defect states are exponentially 
localized and well separated from the rest of the energy spectrnm. Each defect can 
therefore be considered as a qnantnm dot (QD). 

Two snch localized non-degenerate defect states can be thonght of as a basic nnit of 
spin-based information processing. In |9], it was shown that a system of two vacancies 
can be described as a donble qnantnm dot (DQD) system. Fnrther, it was shown that 
the exchange conpling in snch systems can be tnned nsing a magnetic held. 

Snch in situ tnning of the magnetic properties is essential for the realization of 
spintronic devices. In this work, we propose a complementary approach to achieve 
snch tnning based on strain. The electronic properties of graphene are sensitive to 
strain [T^ l2n| [2T], and it is known that GNRs have a strain-tnnable band-gap |22[ |2T]. 
However, there are only a few stndies on the control of magnetism by strain, 
which so far focns on graphene qnantnm dots|23l jH], adatoms |25l |HS] and magnetic 
impnrities|271|28]. The inhnence of compressional strain on vacancy-indnced magnetism 
in bnlk graphene was stndied in Ref. |29j . 

We nse AGNRs as a model system, and characterize the effect of strain on single 
and donble vacancies in snch systems. For a single vacancy we calcnlate the exchange 
splitting within a mean-held approach to the respective Hnbbard model and show that 
it follows the strain-dependence of the band-gap. In the case of two vacancies, we nse 
the localized defect states to setnp a donble qnantnm dot (DQD) modeip, ED]- From 
this model we calcnlate the singlet-triplet energy diherence and hnd that it depends 
exponentially on the prodnct of the distance of vacancies and the applied strain. Onr 
resnlts snggest that strain can be nsed to non-invasively manipnlate the properties of 
magnetic point defects in graphene. This is expected to be of crncial importance for 
spintronic applications. In particnlar, for applications consisting of arrays of snch DQDs, 
strain provides a possibility to access and tnne each DQD individnally. Moreover, by 
combining vacancy-indnced magnetism with mechanical resonators, magneto-mechanic 
devices based on this principle are within reach. 
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We consider ribbons with armchair edges as sketched in Fig. In general, the 
defects can be either (single-atom) vacancies or chemisorpted hydrogen atoms jTT], which 
both effectively lead to the removal of Pz orbitals. In the following, “vacancy” will 
refer to both types of defects. We note that removing a carbon atom also leads to 
an nnpaired cr-electron. Experimental resnlts[3T| snggest that this electron remains 
nnbonded withont hybridizing with the extra 7r-electron. Thns, both electrons provide 
independent contribntions to the total magnetization. As in previons stndiesjlSl [12] we 
only consider the 7r-electrons in the following. 

Uniform strain is applied in the direction parallel to the armchair edges. We inclnde 
the Poisson effect (narrowing of the ribbon npon stretching), and write the elements of 
the strain tensor m = [uxx,Uxy',Uxy,Uyy] as Uxx = u, Uyy = —vu and Uxy = 0, where v 
is the Poisson ratio u the applied strain. The bond-vectors Vik connecting atom i and 
atom k are transformed as > (I -|- m) ■ with I being the 2x2 identity matrix. 
For a realistic interatomic potential, this relation is slightly modified|32|. and the edges 
of the vacancies will be reconstrncted |TT] . These effects depend on the details of the 
interatomic potential, and are not inclnded here. 

To model the electronic properties of the nanoribbon we use a Hubbard model, for 
which the total Hamiltonian is H = + Here, cj^ creates an electron at 

site i with spin s = {t, and U denotes the interaction strength. The non-interacting 
part, is given by the usual tight-binding model with nearest neighbor hopping 
(see for instance [33]). Expanding the hopping amplitude to{r) and the electron-ion 
potential Vo(r) up to first order in small displacements of the atoms, yields the following 
Hamiltonian |34[ l22] l35] 
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Figure 1. Structure of an armchair nanoribbon with two vacancies in a “tail to tail” 
configuration. The missing atoms belong to different sub-lattices and are indicated by 
the blue crosses. The ribbons consist of W x W cells with W counting the cells in 
x-direction and W in y-direction. In the lateral direction periodic boundary conditions 
are imposed at the edges. 
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Figure 2. Strain dependence of (a) the (average) defect energy (edef.,t +£def.,i)/2 
of the localized state, (b) the normalized participation ratio rj{u)/g{Q), (c) the band 
gap Eg and (d) the exchange splitting edef.,t ~ ^def.,^ for AGNRs of different widths. 
Inset in (b) shows the dependence of g{0) on the ribbon width. Symbols denote results 
obtained from a TB calculation and dashed lines indicate behavior according to Eq. 
([|). Here Ni = 100 (for [/ = 0), 7Vi = 60 (for U = IMq) and N 2 = 7,13,19. 


The first sum runs over all atoms and the other sums are restricted to nearest 
neighbors of the atom at site i. The on-diagonal contribution is given by the deformation 
potential and its strength g/rQ = drVo\rQ. The modification of the hopping due to 
displacement of the atoms is determined by fd/ro = —c?r^o/^o|ro- fhe following we use 
to = 2.8 eV, U = 1.6to j36], c/ = 4 eV, /3 = 3.37 |20| and i/ = 0.2|3Z|. 

3. Results 

3.1. Single vacancy. 

First we consider the case of a single vacancy at the center of the ribbon for U = 
0. Since the corresponding localized defect state is always located at the Dirac 
point [13], we expect its energy to follow the behavior of the deformation potential, 
ud = ^ {uxx + Uyy). From Fig. j^a) we infer that this is true for strains up to ~ 3%. 
The localization of the defect state is quantified by the inverse participation ratio 
3 = Hi where fji are the components of the normalized eigenfunction of the defect 
state. For a maximally localized state p = 1, while r/ —)■ 0 for completely delocalized 
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states and large systems. In Fig. |^b) the dependence of rj on the applied strain and on 
ribbon width (inset) is shown. The localization at zero strain decreases with increasing 
width|T5]. For a given width, rj shows a non-monotonic strain-dependence, with a 
minimum that shifts toward higher strains for decreasing ribbon widths. This behavior 
can be understood in a qualitative way by recalling that the electronic energy in the 
Dirac picture is given by = hvpk, where up is the Fermi velocity and k is the electron 
momentum which has the unit of inverse length[22j. The relevant energy scale is set 
by the band-gap Eg, and the corresponding length scale is consequently set by hvp/Eg. 
This suggests that the localization length scales inversely with the band-gap. For a 
semiconducting AGNR the band gap is approximately given hym 


En ~ Sfn min 

^ n=0,l,... 


TT 




in-1/3) 


2 


(1 4- n)u , 


( 2 ) 


which yields the characteristic zig-zag behavior shown as dashed lines in Fig. [^c). Note 
that for the chosen widths the gap is determined by n = 0 (for N 2 = 19 the next 
sub-band starts to contribute around u = 0.04). So indeed, comparing the behavior of 
the participation ratio with the dependence of the band gap on the strain, we find that 
ri{u)/ri{f)) Eg{u)/Egif)), shown as dashed lines in Fig. |^b). 

In the interacting case, Lieb’s theoremj38] predicts that the ground state has total 
spin 1/2. The energies £def.,s are, in general, no longer degenerate and the exchange 
splitting Aex. = £def.,t ~ ^def., 4 . is finite. In the ground state only one of the spin states is 
occupied, which results in a finite magnetization. We have verified this by performing 
mean field calculations for the Hubbard Hamiltonian H. In Figs. |^a)-(c) we show results 
for the average defect energy (£def.,t + ^def., 4 .)/ 2 , the inverse participation ratio and the 
band gap. One can see that the observations made for U = 0 also hold in the interacting 
case. In particular, the relation t]{u)/t]{0) ~ Eg{u)/Eg{0) is still valid. Figure |^d) 
displays the strain dependence of the exchange splitting. Up to a certain strain it 
monotonically decreases and then increases again. Comparing with the band gap shows 
that also Agx. follows the strain dependence of Eg given by Eq. ([^. Altogether, we have 
shown that strain can be used to modify the (localization) properties of the defect state 
and to tune the exchange splitting. 


3.2. Two vacancies. 

Having characterized the influence of strain on the magnetic properties of a single 
vacancy, we now add a second vacancy such that their center of mass is at the center 
of the ribbon. This construct can be thought of as a simple realization of a spin qubit. 
We only consider vacancies with missing atoms in different sub-lattices and in the “tail 
to tail” configuration (see Fig. [^. We hrst characterize the double-vacancy system in 
the non-interacting case, and use the obtained results to construct a double quantum 
dot model for the interacting system. 

In the non-interacting case, we hnd two defect states close to the band center [15]. 
Linear combinations of the corresponding eigenstates yield states localized at the left 
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Figure 3. Strain dependence of (a) the average energy Sc and (b) the energy difference 
of the two defect states, and (c) the normalized participation ratio ri{u)/ri{0) for 
an AGNRs with two vacancies with distance d. Symbols denote results obtained 
from a TB calculation with Ni = 100 and N 2 = 7. The dashed lines in (a) 
and (c) indicate behavior according to itD = 3g/2(l — v)u and Eg{u)/Eg{0) with 
Eg from Eq. ([^, respectively. The dashed-dotted lines in (b) show Ae{u) = 
Ae(0) exp + v)ud^ and the inset in (c) displays 77 ( 0 ) vs the distance d. 

and right vacancy, respectively. Figs. |^a) andj^b) show the average energy of the two 
states and their energy difference Ae: for a ribbon with N 2 = 7. The average energy again 
approximately follows the behavior of the deformation potential and increases linearly 
with the applied strain. The energy difference at zero strain decreases exponentially 
with increasing distance of the vacancies [15]. For a fixed width, the energy difference 
increases monotonically with increasing strain. We hnd that it is approximately given 
by Ae{u) = A£(0) exp , where c ~ y/3. However, as we have seen in 

Fig. He ) the band-gap decreases with strain and therefore at some point the defect 
states start to overlap with the continuum. For the strains shown in Fig. |^b) this 
only happens for the smallest distance (shown as blue circles), since here the splitting 
is largest. The inverse participation ratio, shown in Fig. |^c), displays a pronounced 
decrease when the defect states are shifted into the continuum (blue circles). Otherwise 
T] is found to depend only weakly on the distance (see inset of Fig. [^c)) and it decreases 
approximately linearly with the applied strain, but with a distance dependent slope. For 
vacancies which are far apart the behavior is well described by the ratio Eg{u)/Eg{0), 
as expected for nearly independent vacancies. 

3.3. Double quantum dot model. 

According to Lieb’s theorem|3S] the ground state of the ribbon for 17 > 0 with two single¬ 
atom vacancies in different sub-lattices has total spin zero 5 = 0 (singlet state) |15j. In 
this case the highest occupied state, which corresponds to the lower lying defect state 
is filled with two electrons of opposite spin. To estimate the energy gap to the 
triplet states (5 = 1), we only consider the defect states and take them as highest 
occupied and lowest unoccupied orbitalsjSni Ej- By forming linear combinations of 
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we obtain the states localized at the respective vacancy. Dehning corresponding 

annihilation operators Ons = where n = L,R, we can write the Hamiltonian 

for the double quantum dot model as 

-f^DQD = X! + ^R®Rs®Rs “ t + OrsCILs)) + f4ff X! • (3) 


Here, = ^r equals the average frequency of the two defect states, t = Ae/2 
quantihes the tunnel coupling of the local states and f/gfr = Urj is the local interaction 
strength. 

In the local basis |sl,sr), where sl/r = denotes the spin state of the 

left and right vacancy, the singlet state of the two-state model is given by \I' 5 =o = 
(It,!) - \iA))/\/^ and the triplet states are + I4t))/^, ^s=i = |t,t) 


and 'hgJi = li,i). The respective energies are found from Eq. (3) to be 


Es=o — 


t/eff - + 16t2 


= J, 


pi{o) _pd) _p(~i)_fi 

Es=i — Es=i — — U . 


(4a) 

(4b) 


For |f| f/gff the exchange energy can be approximated by J —4t^/f/eff- 

Using the data shown in Fig. [^we calculate the exchange energy J as a function of 
the vacancy distance and the applied strain. The resulting behavior is shown in Fig. 
Note that we plot J vs du. For zero strain \J\ decreases exponentially with increasing 
distance. Moreover, we hnd that \J\ is exponentially increasing with increasing du, which 
is consistent with our observation that Ae oc exp{2Kdu) with k + /2rQ. Since 

the participation ratio depends only weakly on the vacancy distance, but the energy 
difference strongly decreases, the behavior of J approaches the strong interaction limit 
|t| Ueff for increasing distance. Thus the exchange interaction can be tuned over a 
wide range by adjusting the distance and/or by applying strain to the ribbon. 


4. Conclusions 

In summary, we have characterized the influence of geometry and strain on the electronic 
and magnetic properties of semiconducting AGNRs with one and two vacancies. For a 
single vacancy we hnd that the degree of localization and the spin-exchange splitting 
follows the strain-dependence of the band gap, which is non-monotonously changing. 
This simple relation connecting the intrinsic vacancy induced magnetism to strain via 
the band-gap has, to our knowledge, not been reported previously. It implies that the 
spin-exchange can be changed to large extent by applying strain. Further, the singlet- 
triplet splitting in the two vacancy system depends exponentially on the product of strain 
and inter-vacancy distance. This sensitivity with respect to strain may be exploited in 
the development of quantum-information and sensing applications. In an array of such 
vacancy-induced DQDs, the strain-tuning of the exchange interaction can in principle be 
used to gain individual control of the magnetic properties. Combining strain-sensitive 
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Figure 4. Strain and distance dependence of the exchange interaction J for an 
AGNR {Ni = 100 and N 2 = 7) with two vacancies with distance d. Symbols denote 
results obtained from the DQD model Q and data shown in Fig. Full lines show 
the strong interaction approximation J « —At'^/Uee and dashed-dotted lines denote 
J{u) = J(0) exp ^2-\/3^(l -f v)ud^ . 


magnetic defects with mechanical resonators provides a viable route toward magneto¬ 
mechanical devices|39] and might make it possible to couple vacancy-induced DQDs 
over long distances^]]. 
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